Pacman percolation: a model for enzyme gel degradation 
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We study a model for the gel degradation by an enzyme, where the gel is schematized as a cubic 
lattice, and the enzyme as a random walker, that cuts the bonds over which it passes. The model 
undergoes a (reverse) percolation transition, which for low density of enzymes falls in a universality 
class different from random percolation. In particular we have measured a gel fraction critical 
exponent /3 = 1.0 ± 0.1, in excellent agreement with experiments made on the real system. 

PACS numbers: 64.60.Ak, 82.70.Gg, 87.10.+C 



The extracellular matrix (ECM) is a gel composed by 
various proteins, including collagen, elastin, fibronectin 
and laminin, connected to form an elastic network that 
extends macroscopically. This gel is normally imperme- 
able to cell passage, and ensures organ integrity by insu- 
lating organs and preventing cell dissemination. More- 
over it is the support of cell adhesion and regulates cell 
proliferation, differentiation and locomotion. During spe- 
cific processes, the ECM can be degraded by a variety of 
proteolytic enzymes, especially metalloproteinases, that 
catalyze the hydrolysis of the peptide bonds, increas- 
ing the permeability of the ECM to the passage of cells. 
This degradation process can at some point solubilize the 
gel, realizing a reverse "gel-sol" transition, and bringing 
the ECM to a liquid state, in which cells are no longer 
confined and can freely diffuse. This solubilizing transi- 
tion is especially connected with tumor invasiveness, in 
which some cells access the lymphatic and blood circula- 
tion, and disseminate to distant organs (metastasis). In 
this view, beyond the biochemical processes involved at 
molecular level, the understanding of the physical mech- 
anisms of the ECM degradation is of great importance. 

The passage from a liquid to a gel state, is a critical 
phenomenon 0, 0] j in which soluble monomers bind to 
form larger and larger clusters. At some point, when the 
bond density p becomes grater than some critical thresh- 
old p c , an "infinite" cluster extending macroscopically 
is formed, and the network becomes a solid gel with an 
elastic response. The reverse transition, in which bonds 
are removed, and the system goes from a gel to a liquid 
state, can be clearly described in the same framework. 
The theory of critical phenomena predicts that, near the 
transition, the macroscopic quantities describing the sys- 
tem are related to the distance from the transition (p—p c ) 
by power laws. The average cluster diameter diverges as 
\p — p c \~ u , while the weight average mass as \p — p c | -7 - 
The viscosity diverges as \p — p c \~ k below the transition 
(p < p c ), and stays infinite above it, while the gel frac- 
tion (the density of the infinite cluster) and the elastic 
modulus, that are zero below the transition, grow above 
it respectively as \p — p c \@ and \p — p c \ l - It is important 
to point out that exponents v, (3, 7, etc. . . are universal, 
that is they do not depend on the microscopic details of 
the system, but only on characteristics like the dimen- 



sionality of the system, or whether or not there is a long 
range correlation in the distribution of the bonds. 

Recently, a series of interesting experiments have been 
realized to study the in vitro degradation of protein 
gels by exogenous proteinases, under cell-free conditions 
0.0 ■ In particular, in 3] it was shown that a gel-sol 
transition adequately describes the degradation of the 
gel. Two kinds of gel, fibronectin and ECM gel, and three 
kinds of enzyme, thermolysin, trypsin and proteinase K, 
were used in various combinations. An enzyme solution 
was added to a certain quantity of gel, and the solu- 
bilized fraction X so \{t) of peptides was measured as a 
function of time. The balance between enzyme diffusion 
inside the matrix, and proteolysis reaction, was such that 
the enzyme reached an homogeneous distribution inside 
the gel before the reaction had proceeded to a notable 
extent, so that the experiments studied the volume, and 
not surface, degradation of the gel. The gel-sol transition 
is reached at some time t c , when X so \(t) becomes equal 
to one, or when the gel fraction X gc \(t) = 1 — X so \(t) 
becomes zero. With various gel-enzyme combinations, 
and different enzyme concentrations, it was found that 
A gc i oc \t - t c f, with ~ 1, for t < t c . For ECM gel 
and trypsin for example (3 = 1.01 ± 0.03 3]. It may be 
assumed that the density of hydrolyzed bonds is a linear 
function of time,(p — p c ) cx (t — t c ), at least near the tran- 
sition threshold, so that A gc i oc \p — pd 13 , with (3 ~ 1, for 

P > Pc- 

This result is quite unexpected, because sol-gel tran- 
sition is usually well described by random percola- 
tion, which is obtained when each bond between two 
monomers is present with probability p, and there is no 
correlation between different bonds. Random percolation 
in three dimensions gives a critical exponent f3 = 0.41, 
very different from the one measured in the gel degra- 
dation experiments. In Ref. 3] different possible expla- 
nations of this discrepancy are proposed, including the 
fact that the proteins of the gel may be considered as in- 
terpenetrated polymer coils, that not all the solubilizcd 
fraction may be measured in the experiments, or that 
bonds of the gel phase may be catalyzed more efficiently 
than those of the liquid phase. 

Another possibility, in order to explain the change in 
the universality class with respect to random percolation, 



2 




bond-hydrolyzing enzyme has done 12000 steps, and the den- 
sity of non-hydrolyzed bonds is p — 0.42. Bonds are spatially 
correlated also at large distances, b) The two-dimensional 
random percolation model with the same bond density. 

is the presence of a long range correlation in the distribu- 
tion of non-hydrolyzed bonds g|. The correlation func- 
tion is defined as G(|r|) = (p(r')p(r'+r)) — (p(r')) 2 , where 
p(r) is the density of bonds, and the average (•) is done 
over the reference position r'. It was shown by Weinrib 
and Halperin [6( that if the correlation obeys a power 
law G(r) ~ r~ a at long distances with a < d (where d 
is the dimensionality of the system), then the percola- 
tion transition falls in a universality class different from 
the random percolation, in particular with a correlation 
length exponent v = 2/a. In the case of enzyme catalyzed 
gel proteolysis, the density of non-hydrolyzed bonds is 
certainly correlated, and if the enzyme concentration is 
small it may be correlated also at large distances. 

In this paper, we study a very simple model, which we 
call " pacman percolation model" , in which the protein gel 
is schematized as a cubic lattice oi N — L 3 sites, where 
each site represent an (exavalent) monomer. At time 
t — all the bonds between nearest neighbor monomers 
are present. One enzyme is then introduced in the sys- 
tem, that at each step goes from one site to a nearest 
neighbor site, chosen randomly between the six possi- 
ble neighbors, and hydrolyzes (deletes) the correspond- 
ing bond if not yet hydrolyzed Periodic boundary 





pacman percolation 
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random percolation 


Pc 


0.139 ±0.001 




0.2488 
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1.8 ±0.1 




0.88 





1.0 ±0.1 


1.0 ±0.1 


0.41 
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3.4 ±0.2 




1.80 


/' 


3.5 ±0.1 




2.0 


s 


1.1 ±0.1 




0.73 



TABLE I: Percolation density and critical indices in the pac- 
man percolation model, and in random percolation, in three 
dimensions. 



conditions are chosen. In Fig. QJi it is shown the two- 
dimensional version of the model, after the enzyme has 
walked around for some time (roughly at the percola- 
tion threshold). Note how the remaining non-hydrolyzed 
bonds are spatially correlated, with respect to a random 
percolation model (Fig. U>). At each time step, there 
will be a distribution of clusters, where two sites belong 
to the same cluster if there is a path of non-hydrolyzed 
bond between them. We measure then, as a function of 
the density p of bonds: a boolean variable equal to one if 
there is a percolating cluster, to zero otherwise; the size 
of the percolating cluster, if any; the mean cluster size, 
that is n s s 2 , where n s is the number of clusters of 

size s, and the percolating cluster is excluded from the 
sum. We perform the experiment many times with differ- 
ent random numbers, and average the above mentioned 
quantities, obtaining respectively the percolating proba- 
bility II(p, L), the percolating cluster density p{p, L), the 
mean cluster size x(p, L), where we have explicited the 
dependence on the size L of the lattice. 

From these quantities, it is possible to evaluate the per- 
colation density p c and the critical exponents z/, j3 and 
7 0. Plotting the percolating probability II(p, L) as a 
function of p for different lattice sizes L, it is possible to 
measure the percolation threshold density p c as the point 
in which the different curves intersect, for L — » oo. In 
Fig- Glit is shown the measured H(p, L) for cubic lattices 
of size L = 30, 40, 50, and 60. The intersection point (up- 
per inset) is p c = 0.139±0.001. Plotting then U(p 7 L) as a 
function of (p — p c )L 1 / lJ ', one can measure the correlation 
length exponent v as the value that gives the best col- 
lapse of the curves. In the lower inset of Fig.|5|it is shown 
the obtained data collapse, that gives v = 1.8 ± 0.1. The 
confidence interval is defined by looking when the curves, 
in the interval shown, do not collapse anymore within the 
error bars. In Fig. [3| it is shown the percolating cluster 
density for the same system sizes. Plotting p(p, V)L^I V as 
a function of (p — p^L 1 '", one can measure the exponent 
f3 (see inset), We find /3 = 1.0 ± 0.1, in excellent agree- 
ment with the experimental result of Ref. Q . Finally, in 
Fig.^we measure the mean cluster size exponent 7, find- 
ing 7 = 3.5 ± 0.2. Note that v, (3 and 7 satisfy well the 
hyperscaling relation 2/3 ± 7 = vd, expected on general 
grounds In Tab. [I] the critical exponents found are 
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FIG. 2: Percolation probability II(p, L) as a function of the 
bond density p, for cubic lattices of size L = 30, 40, 50, 
60. Upper inset: the point of intersection of the curves. 
Lower inset: data collapse obtained plotting II(p, L) versus 
(p-p c )L 1/v , with v = 1.8. 
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FIG. 3: Density p(p, L) of the percolating cluster as a function 
of the bond density p, for the same lattice sizes of Fig. [5] 
Inset: data collapse obtained plotting p{p,L)L l3 ^ v versus (p — 
p c )L 1/v , with v = 1.8 and = 1.0. 



compared with those of random percolation. These re- 
sults show that the model of "pacman percolation" falls 
into another universality class with respect to random 
percolation. 

We have also tried to verify the relation predicted by 
Weinrib and Halperin |6| between the exponent v and 
the power law governing the decay of correlations. In 
Fig. El the correlation G(\i — j\) — (rurij) — (ni)(rij) in 
the occupation of the bonds i and j, where \i — j\ is 
the distance between the centers of the bonds, is shown 
for a system of size 100 3 at the percolation threshold 
p = 0.139. The correlation obeys a power law G(r) ~ 
r~ a with a = 1.15±0.05, with an exponential cut-off, 
presumably due to finite size effects, at distances larger 



FIG. 4: Mean cluster size x(Pi L) as a function of the bond 
density p, for the same lattice sizes of Fig. [5] Inset: data 
collapse obtained plotting x{p,L)L ^ v versus (p — p^L 1 ^, 
with v = 1.8 and f3 = 3.4. 
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FIG. 5: Spatial correlation G(r) in the occupation of the 
bonds, for a cubic lattice of size L = 100. 

than r ~ 30. The relation v = 2/ a predicted by Weinrib 
and Halperin, is quite well verified within the errors. It 
has been recently argued Q that for such a model the 
correlation between bonds should decay as 1/r, implying 
a = 1 and v = 2. The prediction, however, is valid only if 
some conditions are verified, such as long times and large 
distances. The discrepancy between this prediction and 
our results may be due to the fact that these asymptotic 
regimes are not reached in our simulations. 

To complete our study, we have analyzed the criti- 
cal indices of the conductivity in the random resistor 
and conductor-superconductor networks. In the first case 
each present bond of the model is substituted with a re- 
sistor of unitary conductance, while absent bonds have 
zero conductance. The total conductivity £ of the model 
is then measured as a function of bond density, and it is 
zero for p < p c , while it grows as \p — p c | p for p > p c . 
Using finite size scaling as usual (see Fig. @) we find 



4 




0.1 0.125 0.15 0.175 0.2 



P 

FIG. 6: Conductivity of the random-resistor network. 
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FIG. 7: Conductivity of the conductor-superconductor net- 
work. 

fi = 3.5 ± 0.1. In the second case each present bond 
of the model is substituted with a superconductor of in- 
finite conductance, while absent bonds are substituted 
with resistors of unitary conductance. In this case the 



total conductivity E diverges as \p — p c \ s for p < p c , and 
stays infinite for p > p c . In this case we find s = 1.1 ±0.1 
(see Fig. 0) 

It was proposed some time ago 0, |^ , that the expo- 
nents p, of the resistor network, and s of the conductor- 
superconductor network, could be in correspondence re- 
spectively with the exponents t of the elastic modu- 
lus, and k of the viscosity. The first correspondence 
is based on the analogy between the equilibrium con- 
dition of sites under elastic tension (JA fcij( r i — r j) = 0) 
and the Kirchhoff equations for the current conservation 
(zCj a ij{Vi = 0)- However, the electrical problem 

is a pure scalar problem, while the elastic one has a vec- 
torial nature, so the macroscopic elastic response is in 
principle determined by more than one elastic modulus. 

It is worth to notice that, if the density of enzyme in- 
creases, the process of degradation of the gel becomes 
closer to a random degradation. Indeed, in the limit in 
which each bond is hydrolyzed by a different enzyme, 
there will be no correlation between them. This is con- 
firmed by our simulations. We have analyzed the behav- 
ior of the model with a finite density pE of enzymes (de- 
fined as the ratio between the number of enzymes present 
and the number of sites of the lattice), and we have 
found that for densities greater than pe — 0.8 the system 
falls in the universality class of random percolation. For 
p E = 0.8 we find v = 0.88 ± 0.02 and (3 = 0.41 ± 0.01, in 
excellent agreement with random percolation 0. 

In conclusion, we have used a percolation model to 
study the degradation process of ECM due to the ac- 
tion of enzymes. Our results shows that, for low density 
of enzymes, our model belongs to a different universal- 
ity class from random percolation. The change in the 
critical indices may be due to long range correlation. If 
the density of enzymes is sufficiently high, the correla- 
tion between bonds disappears and there is a crossover 
to random percolation. 
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